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Berry Curvature of the Dirac Particle in a-(BEDT-TTF)2l3 
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We examine several properties of the Berry curvature for the organic conductor a-(BEDT- 
TTF)2l3 consisting of four bands, which exhibits a zero-gap state with Dirac cones. By adding 
a small potential acting on two molecular sites, which breaks the inversion symmetry, it is 
shown that the curvature for the Dirac particles displays a pair of peaks with opposite signs 
and that each peak increases with decreasing potential. The Berry curvature originating 
from the property of the wave function is analyzed using a reduced Hamiltonian with a 
2x2 matrix based on the Luttinger-Kohn representation, which describes a pair of Dirac 
particles between the conduction band and the valence band. Two types of velocity fields in 
the reduced Hamiltonian, whose vector product gives the Berry curvature, rotate around the 
Dirac point as a vortex. It is also shown that the other bands exhibit another pair of peaks 
of Dirac particles with a tendency toward merging. 

KEYWORDS: Berry curvature, Dirac particles, Dirac point, zero gap, Luttinger-Kohn repre- 
sentation, a-(BEDT-TTF)2l3, organic conductor 

1. Introduction 

The metal-insulator transition in the two-dimensional organic conductor Q!-(BEDT- 
TTF)2l3^) (BEDT-TTF=bis(ethylenedithio)tetrathiafulvalene) has been studied extensively 
since the discovery of electric conductivity, which exhibits a noticeable dependence on both 
pressure and mag netic field.2) With increasing uniaxial pressure Pa along the stacking axis, 
the insulating state changes into the superconducting state, and is followed by a narrow- 
gap semiconducting state. ^-^ The Hall conductivity at Pa = 10 kbar decreases rapidly with 
decreasing temperature.'^^ Using the transfer energy obtained from the structure analysis, 
the successive transition has been studied theoretically in terms of an extended Hubbard 
model. ^) In addition to the vanishment of the density of states in the narrow-gap state, it 
has been demonstrated that the zero-gap state exhibits a contact point followed by a pair 
of Dirac cones. The zero-gap state was confirmed by first-principles calculation.^' Such a 
Dirac particle, can explain the anomalous temperature dependence of resistivity^*^) and Hall 
coefficient.^^) 

In Fig. 1, the crystal structure of a-(BEDT-TTF)2l3 is shown where there are four sites. 
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Fig. 1. (Color online) Crystal structure on two-dimensional plane with four molecules A(l), A'(2), 
B(3), and C(4), in the unit cell where the respective bonds represent seven transfer energies 



i.e., A, A', B, and C, in the unit cell with an inversion symmetry between the A and A' sites. 
For the insulating state, which is found at ambient pressure, the symmetry is broken since 
the electron density at the A site is different from that at the A' site owing to a stripe charge 
ordering perpendicular to the stacking axis.^^"^^^ At high pressures where the zero-gap state 
emerges, the symmetry is restored, and the electron density at the A site becomes equal to 
that at the A' site.^*^) The zero-gap state is robust to the variation in pressure, although the 
contact point moves in the Brillouin zone (BZ).^) Small amounts of anion potential which 
are the same for the A and A' sites but are different for the B and C sites, do not break 
the contact point of the Dirac cone.^^^ These site potentials for B and C also stabilize the 
zero-gap state by reducing the overlap between the conduction and valence bands. A pair 
of Dirac cones located at the incommensurate wave vector in the BZ is expected to vanish at 
high pressures by merging at a symmetric point (e.g., T point). ^^"^^^ Furthermore, it has been 
shown that the four sites, A, A', B, and C provide the respective role for the tilted Dirac cone, 
and that each site gives a difference in the local density of states, as shown in the temperature 
dependence of the Knight shift and Tf ^ of 

Such a Dirac particle can be represented in terms of a reduced Hamiltonian with a 2x2 
matrix instead of the original 4x4 matrix Hamiltonian, where the reduced Hamiltonian based 
on the Luttinger-Kohn (LK) representation,^^) gives the tilted Weyl equation^^'^^'^^^ for the 
Dirac particle. The reduced model is successful to examine the novel property of transport 
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phenomena/^' ^^'■^^^ which exhibit several differences from those of graphene,^^'^^^ owing to 
tilting. 

Although much progress has been achieved on the basis of the zero-gap state ^29-31) ^j^g 
condition for the Dirac particle is not yet clear since a contact point is formed owing to 
the accidental degeneracy in the energy band.^^^ In the present study, as the next step in 
understanding the Dirac particle, we examine the Berry curvature'^^"^^-* for the a-(BEDT- 
TTF)2l3 salt, which comes from the property of the wave function. The curvature is also a 
useful tool for searching for the particle in other organic conductors. The Dirac particle has 
been maintained using a fact that the gap between the first band and the second band is zero 
within numerically accuracy. However, this is not always the case for the Dirac particle since 
there is also a Dirac particle with a gap, e.g., a massive Dirac particle in the low-energy model 
of boron nitride, and also in the charge ordered state of a-(BEDT-TTF)2l3 at pressures 
slightly lower than those of the zero-gap state. 

In §2, the formulation for the Berry curvature of a-(BEDT-TTF)2l3 is given using a 4x4 
Hamiltonian. To calculate the Berry curvature with a finite magnitude around the contact 
point, we add a small potential that differs between the A and A' sites, based on the previous 
work.^^) Furthermore, a reduced Hamiltonian of a 2x2 matrix is examined by calculating the 
analytical property of the curvature. In §3, the Berry curvatures of the Dirac particle with 
multibands are studied numerically by taking a typical pressure. Pa =6 kbar. Velocity fields 
and Berry curvature are also studied. The summary and discussion are given in §4. 

2. Formulation and Hamiltonian 

2. 1 Berry curvature of multibands 

We consider a two-dimensional electron system consisting of M molecular sites in a unit 
cell where M = 4 in the present case. The Schrodinger equation is given by 

H{k)\n{k) >= En\n{k) > , (1) 

where k = {kx,ky) is a two-dimensional wave vector, and n(= 1,2,- •• , M) denotes a band 
index. The quantities H{k), En{k), and \n{k) > are the Hamiltonian with an MxM matrix for 
the M sites, eigenvalue (band energy), and eigenfunction (wave function), respectively. When 
k{t) varies slowly with time t, the equation for motion is written as 

ihdt\n{k{t)) >= H{k{t))\n{k{t)) > . (2) 
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When k moves along a closed loop C, the state \n[k) > after returning to the original location 
is replaced by e*'^"|n(fc) >. The quantity 7„ is the Berry phase given by^^^ 

7„ =i / {n{k) • Vfcn(fc)) dk= [ Bn{k) ■ dS, (3) 
Jc Js 

where S is a region enclosed by the loop C. The vector field Bn{k) is the Berry curvature 
expressed as 

Bn{k) = -Im (Vfcnl x | V^n) , (4) 

which is oriented perpendicularly to the kx-ky plane. Noted that -B„ = for a single component 
owing to < n\n >= 1. 

The Berry curvature can be rewritten as^'^^ 

B„(fc) = -In, <"l^-y_-<"'l7>^l"> . (5) 

which is used for the numerical calculation. Note that 

4 

n=\ 

for an arbitrary k since (?7-|Vfc-ff|m) is complex conjugate to (m|Vfci7|n), and the double 
summations of n and m result in the vanishment of the imaginary part. This suggests a possible 
pair of Dirac cones located between two neighboring energy bands, because the singularity of 
one band must be compensated with that of the other band. 

2. 2 Hamiltonian for a- (BED T- TTF)2h 

The contact point, which is obtained in a normal state, i.e., without charge ordering, is 
determined essentially by the property of the transfer energy, and the effect of interaction is 
to modify mainly the location of the contact point (see Fig. 5 in ref. 19). Then, by taking 
account of only transfer energy, we examine the Hamiltonian for a-(BEDT-TTF)2l3 given by 

4 N 
a,l3=l o ij 

where a, (3 (=1, 2, 3,4) are indices for the sites of molecules A (1), A'(2), B (3), and C (4) in 
the unit cell, and i,j{= Ri,Rj) are those for the cell forming a square lattice with N sites. The 
quantity aj^^- denotes the creation operator for the electron, and ta^p-ij denotes the transfer 
energy between the neighboring site. As shown in Fig. 1, there are seven transfer energies 
given by tfoi • • • for the direction of the 6(— x)-axis and by tai - ■ ■ tas for the direction of the 
a(?/)-axis. We note that the definition of transfer energies in Fig. 1 (ref. 1) differs from our 
previous one, tpi, ■ ■ ■ ,tcA-'^^ Since the two-dimensional plane of the previous one is reversed 
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compared with that in Fig. 1, the relation is given by tpi = tb2, tp2 = ^bi; ip3 = tb4, tp4 = tbs, 
tci = ta2i tc2 = taSi and tc3(= tc4:) = tai- In Order to calculate the Berry curvature, which 
becomes finite around the Dirac point, we introduce site potentials laa acting on A (1) and 
A' (2) differently in the form of 

4 N 



Q,/3=l i,j 



where 



la 



-L 



2a 



-An 



(9) 



and /30- = lia = 0. Equation (8) is a Hamiltonian of the site potential that creates a gap by 
breaking the inversion symmetry between A and A'. 

Using the Fourier transform aa{k) = J2j (^a,j exp[—ik ■ Rj], the total Hamiltonian is 
expressed as 



site 



H{k) 



J](ai,(fc)t, a2.(fc)t, a3.(fc)t, ai^{k)A Hik) 



a2aik) 
asaik) 
\a4aik) J 



I 



-An 



tcl + tc2G' 

Ao 



(10) 



\tp2 — tpse 



ihx 



ihx 



tp3 — tp2e' 



ihx 



tpie'-^y - fpie*(^'^+''!') 


tr 



r^ihx 



J 



tp2 
pi — tp2e 

tc3 + tc4e~'''y 


(11) 

The lattice constant is taken as unity. By diagonalizing eq. (11), we obtain four bands Ej(k) 
(j = I,-- - ,4), where £'i(k) > ii^2(k) > i?3(k) > £'4(k). For laa = where the inversion 
symmetry between A and A' is maintained, the zero-gap state between i?i(k) and -E'2(k) is 
obtained at =b/co, i.e., with the uniaxial pressure along the a-axis, Pa, being larger than 3 
kbar.'') In the present paper, we introduce a small but finite Aq that induces a small gap 2A 
between two upper bands. We call /cq the Dirac point where fco ^or Aq = gives the contact 
point, and fco for Aq ^ corresponds to the minimum of the gap between Ei[k) and E2{k). 



2.3 Reduced Hamiltonian 

We analytically examine -Bi, which denotes the curvature for the first band Ei. The 
present numerical calculation shows that Bi in eq. (5) is essentially determined by Ei and 
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E2, i.e., the effects of and E4 on Bi are neghgibly small. Thus, eq. (5) may be rewritten 
as 

B(k)=- r im\^kH\2{k)) X (2(fc)|Vfcif|l(fc)) 

{Vi{k) X V2{k)), (12) 



{E2-E^)^ 

where vi and V2 are the velocity fields defined by 

(l(fc)l Vfc//|2(fc)) = vi{k) + iv2{k) . (13) 

Although the velocities vi and V2 depend on the choice of the relative phase between |l(fc) > 
and |2(fc) >, the product of vi x V2 is independent of the phase. 

Prom eq. (13), one can reduce the 4x4 Hamiltonian to the 2x2 Hamiltonian for the energy 
bands Ei and E2. Instead of eq. (12) where the band basis on k is taken, the effective Hamil- 
tonian is expressed using the LK basis, which is fixed at fc^^(= feo)- The general form of the 
2x2 Hamiltonian, H^'^'^{k), is given by 

\fi + if2 -/3 + /0, 

= /O + /iC^l + /20-2 + /3<73 • 

(14) 

The quantity aj denotes the Pauli matrix and the coefficient fj (= fj{k)) is given as the 
function of the two-dimensional wave vector k. Although the term /o gives the effect of the 
tilting for the Dirac cone,^^^ /o is irrelevant to B and is discarded hereafter. Substituting 
eq. (14) into eq. (12), we obtain the Berry curvature as (see Appendix) 

B = (/3(Vfc/i X Vfc/2) + /l(Vfc/2 X Vfc/3) 

+/2(Vfc/3 X Vfc/i)) , (15) 



where E = ^ f( + fi + /|. 

Now, we examine the Berry curvature close to the Dirac point k^ using the LK basis at 
fc^^. The quantity fj{k) is expanded as 

h(k)^vl-q, 

f2{k) ^ vl ■ q, 

h{k) -A, (16) 
where q = k — k^^, v^{= ^kfj) is estimated at fc = k^^, and the linear term in /3 is absent. 



6/18 



J. Phys. Soc. Jpn. Full Paper 

The quantity 2A(< 2Ao), which is the gap between the two bands of Ei and E2, is small 
compared with the bandwidth. From eqs. (15) and (16), the Berry curvature and Berry phase 
are respectively obtained as 

^ = ^ (-? X , (17) 

71= /(e.-B)dg = vr "--;i"i; (18) 

where E = + /I + and denotes the unit vector perpendicular to the qx-Qy plane. 
In deriving eq. (18), the two-dimensional integral is extended to infinity. 

Now, we examine the variation in the velocity fields around the Dirac point for an arbitrary 
fc(^ fcLK-j ^gjj^g (13). By considering the case of 

v^^^{v,Qi), vl^{{),v') and /g ^ -A, (19) 

we obtain (see eqs. (A-9) and (A- 10) in Appendix) 

MOq) = -^-j^=={v^qxex + {v'fqyey) , (20) 

MGq) = -^^==^vv' {-qyEx + q^Sy) , (21) 
„ AW 

B = , (22) 



where E = ^^v^q^ + + A^. These velocity fields rotate around fco as shown explicitly for 
the 4x4 Hamiltonian in the next section. 

Here, we comment on the reduced Hamiltonian of eq. (19) rewritten as 

TTin I ~^ qx-iqy\ 

' = I = qxcri + qycr2 - Aa^ , 

\qx + iqy A / 

(23) 

where v = v' = 1. By introducing a unitary transformation, eq. (23) is transformed as follows. 
The transformation for e"*""^'^/^ with (j) = -k /2 gives 

= e"*'^2</>/2j^(/)gjf72</./2 

qx A - iqy 
A + iqy -q. 



(24) 



and that for q^'^^^I'^ with (f) = tt/2 gives 
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Fig. 2. (Color online) Energy bands of Ei{k) (a) and E2{k) (b) for Pa = 6 kbar and Aq =0.02 eV, 
where the regions for the and ky axes are given as < fc^; < 27r and — tt < ky < n. The arrow 
denotes the location for the Dirac point at kp = (O.STtt, — O.SOtt) where the gap between Ei{k) 
and E2{k) is estimated as 2A ~ 0.02 eV. Two cones are situated symmetrically with respect to 
the X point (tt, 0) and the P point (0, 0). 

(ky kx — 

= kxcri + AfT2 + kyaz ■ 
kx + iA -ky J 

(25) 

For A = 0, eqs. (24) and (25) are reduced to those derived by Katayama et al.^'^) and Kobayashi 
et al./^) respectively, who chose different bases for the L-K state with and V2 being or- 
thogonal. Note that all these cases lead to -B^ = S • = A/(2£^'^) owing to the rotation that 
does not depend on q. 

3. Berry Curvature for Q!-(BEDT-TTF)2l3 

To calculate the Berry phase of the conductor under uniaxial Pa along the stacking axis, 
we use an extrapolation^^ for 1^/3 given by^^ 

tapiPa) = ta(3i0) + C^^Pa . (26) 

The unit of energy is taken as eV and that of pressure Pa is taken as kbar. Equation (26) 
corresponding to Q-(BEDT-TTF)2l3 gives the zero-gap state for Pa > 3 kbar. We perform 
the calculation by taking Pa = 6 kbar and Aq = 0.02 eV. The gap on the Dirac point ko, i.e., 
2A(~ 0.02 eV), is smaller than 2Ao. 

The energy bands Ei and E2 are shown in Figs. 2(a) and 2(b), respectively. The Dirac 
cone exists at k ~ itko with kg = (0.577r, — O.Svr), which is shown by an arrow. A pair of 
cones is seen around £'i(ko) and £'2(ko), and the cones in the same band are symmetric with 
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Fig. 3. (Color online) Berry curvature and the contour plot of —B ■ for the energy band Ei(k) 
on the kx-ky plane corresponding to Fig. 1(a), where denotes the unit vector perpendicular to 
kx-ky plane. The middle of the horizontal axis (vertical axis) in the contour plot denotes the M 
(P)) point. 

respect to the r(=(0,0)) point. For Ei{k), the maximum is seen at the Y(=(0, itvr)) point 
while saddle points are seen for the X(=(7r,0)), M(=(7r, ztyr)), and T points. For E2{k), the 
minimum is seen at the M point while saddle points are seen for the X, Y, and F points. The 
accidental degeneracy, which is found at itfco for Aq — )• 0, is realized as the minimum of Ei(k) 
and the maximum of E2{k). 

The Berry curvature, which is calculated using eq. (5), is shown in Figs. 3(a) and 3(b) 
for the three-dimensional (3D) view on the kx-ky plane (a), and for the contour plot (b). The 
peak around —k^ ( k^) is positive (negative). They are antisymmetric with respect to the F 
point (the middle point on the vertical axis). Since the curvature exhibits a noticeable peak 
close to fco) such a peak may be identified as the Dirac particle instead of calculating the 
contact point. We note that a quantity called a Chern number'^^^ is given by 



owing to the time reversal symmetry, i.e., B{—k) = —B{k), where S denotes the total BZ. 
In Fig. 3(a), the curvature around the Dirac point is slightly anisotropic in the sense that it 
is slightly extended to the y-direction. With increasing pressure, the difference is reduced and 
the curvature becomes almost isotropic at Pa = 10 kbar, where two Dirac particles are well 
isolated. 

Using eq. (13), we examine velocity field as a function of k. In Figs. 4(a) and 4(b), the 
k dependence of the velocities vi and V2 around fco are shown. Since the phases of both 
|l(fc) > and \2{k) > are arbitrary, we choose them so that both of their first components 




(27) 
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Fig. 4. (Color online) Vector fields of velocities — vi (a) and —V2 (b) on the kx-ky plane around the 
Dirac point, feo, (cross) in Fig. 2. 




(b) 



Fig. 5. (Color online) Contour plot of Berry curvature, B{Aq), for energy bands E2{k) (a) and E^ik) 
(b). 



corresponding to the A site are real. We verified that B calculated in terms of vi{k) x V2{k) 
in eq.(12) reproduces well the result in Fig. 3 within numerical accuracy. Both vi and V2, 
which are almost orthogonal, rotate around the Dirac point, indicating the singularity at kQ. 
The vortex structure around the Dirac point is also understood from eqs. (20) and (21). Even 
in the presence of the site potential Aq, which removes the degeneracy of the zero-gap state, 
the existence of such vortex behavior suggests a topological property of the Dirac particle. 

Here, we examine the Berry curvature for other bands of E2{k), £^3(fc), and E4{k). The 
Berry curvature of the second band E2{k) is shown in Fig. 5(a). The peak located at itfco has 
a sign opposite to that of Ei{k). In addition to such a peak, another pair of peaks appears 
close to the P point. The latter one is rather extended to a direction slightly declined toward 
the horizontal axis, suggesting a large anisotropy of the Dirac cone. The anisotropic peak 
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Fig. 6. (Color online) Contour plots of the Dirac eomponent of eq. (29) (a), and the background 
component of eq. (30) (b), for Berry curvature of Ei{k) in Fig. 3. 



disappears for a small < 3 kbar while it becomes rather isotropic for a large Pa- Such a 
behavior resembles the emergence of the Dirac particle in the charge ordered state, which has 
been shown in a- (BEDT-TTF)2l3.^^) In Fig. 5(b), the Berry curvature for E^ik) is shown. A 
pair of peaks close to the F point also shows a sign opposite to that of E2{k). There is another 
peak close to the M point, which compensates for that of i?4(fc) (not shown here). These 
results show that each neighboring band provides a pair of Dirac particles followed by the 
Berry curvature with an opposite sign. When one pair of Dirac particles between neighboring 
bands is found, one may expect another pair between neighboring bands. 

Since the Berry curvature B{k) in Fig. 3 depends on the sign of the site potential =bAo, we 
examine B(k) by dividing it into two components, i.e., -Boirac and -Bbacki which are symmetric 
and antisymmetric, respectively, with respect to Aq. By defining 



B{Ao) = Bi{k,Ao)-e, 
these components are rewritten as 

irac 



B{Ao)-B{-Ao) 



-Bback 



B{Ao) + B{-A, 



(28) 

(29) 
(30) 



The first one is related to the topological property, which does not depend on the details of 
the energy band. The second one depends on the property of the energy band. The curvatures 
given by eqs. (29) and (30) are shown in Figs. 6(a) and 6(b), respectively. The peak height 
of i?Dirac IS 10^ times as much as that of -Bback- When Aq decreases to zero, i?Dirac increases 
rapidly, but -Bback remains almost unchanged. The component -Boirac is located in the narrow 
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Fig. 7. Ao dependence of fiDirac (eq. (31)), f^back, (eq. (32)), fiabs (eq. (33)), and gap 2A . 



region around ±fco) suggesting a character of the Dirac particle. The component -Bbackj which 
is smah and spreads in the entire region of the BZ, is irrelevant to the Dirac particle. Thus, 
the separation of -Boii-ac from B{Ao) gives the intrinsic curvature of the Dirac particle. 

To estimate the magnitude of the peak of B(k, Aq), we examine the following integrated 
quantities of the respective components: 



i^Dirac = — / sgn{ky) 



n 



2tt 

1 



back 



dS sgn{ky 



(31) 



(32) 



2tt Js " 2 

The quantity J^Dirac + ^^back is a Berry phase where 5 in eq. (3) is the region given by — vr < 
kx < vr and < A;^^ < vr. The quantity f^Dirac comes from the contribution of the Dirac cone 
with a singularity, and the second one originates from the property of the conventional band. 
Note that iloirac remains the same but rJback varies when we replace sgn{ky) by sgn{kx) in the 
limit of small Aq. In Fig. 7, l^Dirac and riback are shown as the function of Aq. With decreasing 
Aq, f^Dirac is reduced to 1, since the boundary effect of BZ becomes negligible. This indicates 
the intrinsic property of Dirac cones. We also show the quantity 



1 



|i?(Ao)| + |i?(-Ao) 



Oab. = f dS '"^^"^'"'^^"^"^ (33) 

Note that r^abs is larger than J^Dirac and riback) and becomes even larger than 1 for a small Aq. 
The enhancement of r^abs is attributable to the alternation of the sign of i3(Ao) as a function 
of k, which originates from the background property of the band structure. 
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4. Summary and discussion 

We obtained the following for the Berry curvature of Q-(BEDT-TTF)2l3, which exhibits 
the zero-gap state under the uniaxial pressure. 

(i) In the present study, instead of the zero gap between Ei[k) and E2{k), the Dirac 
particle is examined by calculating the Berry curvature for the 4x4 Hamiltonian. By adding a 
small potential Aq acting on the A and A' sites with opposite signs, which breaks the inversion 
symmetry, the Dirac particle is identified by the pronounced peak of the Berry curvature in 
the same band, and also between neighboring bands. The Berry curvature consists of two 
components, i.e., -Boirac and i?back> where -Boirac (-Bback) is an odd (even) function with respect 
to Aq. The peak in Soirao which is located close to the Dirac cone, is intrinsic to the Dirac 
particle. The quantity -Bback) which is small but is extended in the whole BZ, is irrelevant 
to the Dirac particle. The latter could be the origin of the curvature, which depends on the 
choice of the phase in the transfer integrals of the site Hamiltonian. ^''^ 

(ii) The Berry curvature is calculated in the general case of the reduced 2x2 Hamiltonian. 
Once such Hamiltonian (i.e., the coefficient of Pauli matrices) for other organic conductors 
with possible Dirac particles is obtained on the L-K basis, eq. (15) is applicable to understand 
the characteristic behavior of the Berry curvature. The Berry curvature is determined by a 
vector product of two types of velocity fields, which are given as the coefficients of the Pauli 
matrices, ai and a2, in the reduced Hamiltonian. These velocity fields rotate around the Dirac 
point as a vortex and are determined uniquely, except for the Dirac point. 

(iii) We found a pair of Dirac particles not only between Ei{k) and E2{k) , but also 
between the other neighboring band, although the mutual relation of these Dirac points is 
unclear owing to the accidental degeneracy. In the present case, the curvature between E2{k) 
and E^{k) exhibits a behavior close to merging. It is possible that a pair of Dirac particles 
may exist in many neighboring bands, if a pair of peaks is found in one band. 

Here, we comment on the effect of a short-range repulsive interaction on the Dirac particle. 
In the presence of the interaction that also gives the normal state, i.e., no gap, the critical 
pressure for the zero-gap state {Pa ~ 4 kbar)^^^ is larger than that in the absence of interaction 
{Pa ~ 3 kbar), suggesting that the interaction suppresses the zero-gap state. Such an effect 
appears as the on-site potential acting on the B and C sites. Comparing the Berry curvature 
in Fig. 3(a) with that in the presence of the interaction (not shown here), we found that the 
Berry curvature in the presence of the interaction becomes rather anisotropic and exhibits 
a small tail toward the M-point, indicating a tendency toward merging. Thus, merging may 
occur owing to the interaction, in addition to the variation in the transfer energy.'^) Note 
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that the pair of Dirac particles aheady merges at the M-point in the charge ordering state at 
ambient pressure, which was obtained in our previous study.^^ 

At present, it is not yet clear how Dirac particles emerge in organic conductors, since it 
occurs accidentally in the BZ. However, the calculation of the Berry curvature is a useful tool 
for searching Dirac particles particularly in the multiband system even when the symmetry 
is broken or when the band structure is complicated. Actually, it is shown that the organic 
conductor a-(BEDT-TTF)2NH4Hg(SCN)4 with a complicated band structure exhibits Dirac 
particles with a zero gap under pressure. In the case where the symmetry is broken, the 
formation of massive Dirac particles in the stripe charge ordered state is reported in a separate 
paper. 
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Appendix: Berry curvature and velocity field for general Hamiltonian 

We calculate the Berry curvature for the reduced 2x2 Hamiltonian given by (eq. (14)), 

H^^{k)=( (A-l) 

V/l+^/2 -h ) 

where the quantity aj denotes the Pauli matrix and the energy is measured from the chemical 
potential. The quantity fj (= fj{k)) depends on the two-dimensional wave vector fc(= ko + q). 
The Schrodinger equation for the Hamiltonian (A-l) is written as 

^rcd^^ = ±^^±, (A-2) 
where \I'+ = |l(fc) >o, ^~ = |2(fc) >o- The energy and wave function are calculated as 



E = \/ff + fi + fi, (A-3) 
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The relative phase between ^+ and is undetermined and is chosen such that the second 
component of each state is real and positive. 

Now, we calculate Bi for ^_|_ explicitly using eq. (4). By noting that only the first com- 
ponent contains the imaginary part, the Berry curvature B{k) is calculated as 

Bi = B = -Im fVfc^O* X (Vfc^'+) 



where 



-Im(Vfc/i + iVkh) X (Vfc/i - iVkh) = 2Vfe/i X V^/a , (A-5) 



{2E - /3)/2/i {E - h)fi 

~ 2E\E-h) + V./3} , (A.6) 



and 



- ^2EiE - h) ^- 2EHE - h) ^'f' 

+ " 2i?2(i^_/3)^'=^2 + ^^^Vfc/s}. (A.7) 

Substituting eqs. (A-6) and (A-7) into eq. (A-5), we obtain eq. (15). 

Note that eq. (15) can also be obtained from eq. (5), which is rewritten as 

B = -Im |(*^Vfcif*_) X (vI/t_Vfci7vI/+)-L| = --Lx 

Im{ (*+|cTiVfc/i + CT2Vfc/2 + a3Vfc/3|*-) X (^'-|aiVfc/i + CT2Vfc/2 + (73Vfc/3|^+) }, 

(A-8) 

where (^±| = (/i + if2,g±Y- Equation (A-8) is calculated as eq. (15) by noting 
(^_|cTi|^+) = -(/3/1 - iEf2)/{E.JW^), (^_|a2|^+) = -(/3/2 + iEh)/iE^E^-fi), 



and (^_|a3|^'+) = y^W~fi/E. 

Here, we calculate the velocity field defined by eq. (13). Substituting eq. (A-4) into eq. (13), 
we obtain 

viiOk) =Re{^'lVfcifLK^+} 
-1 



(/3(/lVfe/i + /2Vfc/2) 



-Ul + fh^kh). (A-9) 
V2{ek) =Im{^_Vfcii-LK^^} 
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-1 



-/2Vfc/i + /iVfc/2) , (A-10) 



= } A h^kh + /2Vfc/2 + /aVfc/s) . (A-ll) 

V - fi 

Substituting eqs. (A-9) and (A-10) into eq. (12), one finds that Bi is exactly the same as 
eq.(15). 
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